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Light-assisted reactive collisions between laser-cooled Ba^ ions and Rb atoms were studied in 
an ion- atom hybrid trap. The reaction rate was found to strongly depend on the electronic state 
of the reaction partners with the largest rate constant (7(2) x 10~^^ cm"^ s~^) obtained for the 
excited Ba^(6s)+Rb(5p) reaction channel. Similar to the previously studied Ca^+Rb system, 
charge transfer and radiative association were found to be the dominant reactive processes. The 
generation of molecular ions by radiative association could directly be observed by their sympathetic 
cooling into a Coulomb crystal. Potential energy curves up to the Ba^(6s)+Rb(5p) asymptote 
and reactive-scattering cross sections for the radiative processes were calculated. The theoretical 
rate constant obtained for the lowest reaction channel Ba^(6s)+Rb(5s) is compatible with the 
experimental estimates obtained thus far. 



I. INTRODUCTION 

Following the impressive developments in the generation of translationally cold atoms, molecules and ions [IHS], 
the study of chemical reactions at ultra-low energies has recently emerged as a new exciting research area in chemical 
physics [SHTO]. Over the past few years, ion-atom hybrid traps in which laser pJT[T5] or sympathetically p^BHTS] cooled 
ions in a radio frequency (rf) ion trap are combined with ultracold atoms in a magnetic [12 , optical dipole fi3\ fT8] 
or magneto-optical [iTJ [T4HT7] trap have paved the way for the investigation of sub-Kelvin collisions of atomic and 
since recently also molecular [17] ions with ultracold atoms. At these collision energies, the quantum character of the 
collision can influence elastic \TW as well as inelastic [20l [2Tj scattering cross sections and the validity of ab-initio and 
quantum scattering calculations as well as reactive capture models can be checked accurately [ll| [T4 | \TE\ \T7\ [22l [23l. 
Moreover, the possibility to introduce near-resonant laser fields reveals the prominence of radiation-driven chemical 
processes and highlights the role of electronically excited states in cold collisions through the determination of state- 
specific rate constants [Ml [T5l I22H24] . 

In the present paper, we report on a joint experimental and theoretical study of light-assisted reactive collisions 
between laser-cooled Ba+ ions and Rb atoms extending the data on collisions between the ground-state species 
reported by Schmid et al. [13]. Utilising hybrid-trapping technology average collision energies (^co1i/^b) ^ 600 mK 
were achieved. Rate constants for excited reaction channels populated by near resonant laser fields were determined. 
We find significant enhancement of the rates in electronically excited channels which is rationalised with reference 
to computed potential energy curves and transition dipole moments. Quantum-scattering calculations predict that 
reactive collisions between the ground-state collision partners predominantly lead to the radiative association (RA) 
of molecular ions which is supported by the experimental observation of BaRb+ ion formation. The collision energy 
dependence of the reactive processes is studied and the experimental rate constants are compared with theoretical 
predictions. Finally, the present results on Ba++Rb are compared with and contrasted to the findings in other 
ion- atom hybrid systems. 

II. METHODS 

The experimental setup used for the simultaneous trapping of atomic ions and neutral atoms has already been 
described previously [5l [Ml [23] and only aspects relevant to the present study are described here. 

Ba+ ions were produced by non-resonant photoionisation of a beam of Ba atoms with their natural isotopic 
abundances at the centre of a linear radio frequency (RF) ion trap [5] [H [23]. Trapped i3^Ba+ ions were laser 
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FIG. 1: (a) Laser cooling schemes for ^^^Ba^ and ^^Rb. (b) Superposed false-colour fluorescence images of a Coulomb crystal 
of "«Ba+ ions (blue) and a cloud of ultracold ^^Rb atoms (yellow-red) in the hybrid trap. The asymmetric shape of the ^^^Ba^ 
Coulomb crystal is caused by non-fluorescing sympathetically-cooled Ba^ isotopes, see text for details. 



cooled to secular temperatures of ^ 10 mK using a laser beam at a wavelength around 493 nm slightly red detuned 
from the {6s) ^5'i/2 (6p) transition. A repumper laser at 650 nm was used to repump population on the 

(bd) ^1^3/2 (6p) ^Pi/2 transition (see figure 1(a)). Under these conditions, the Ba+ ions formed Coulomb crystals 
(Fig. [ijb)) [5J which were imaged by collecting the spatially resolved laser-cooling fluorescence of the ions using a 
charge coupled device (CCD) camera coupled to a microscope. 

^^Rb atoms were produced by a getter source and laser cooled and trapped in a magneto-optical trap (MOT) |l4( f23]. 
95% of the MOT laser beam intensity was detuned by ^ 20 MHz from the (5s) ^5'i/2,i^ = 2 ^ (bp) ^P3/2,i^ = 3 
cooling transition of ^^Rb around 780 nm. The remaining 5% were tuned into resonance with the (5^) '^Si/2iF = 1 
{bp) ^P3/2,F = 2 repumper transition (see Fig. [l|a)). The number density {uRb ^ 1 x 10^ cm~^) and temperature 
(T = 90 — 150 /iK) of the atoms in the MOT were determined using standard fluorescence measurement and time- 
of- flight methods, respectively [25l[26]. To avoid multi-photon ionisation of Rb atoms in the excited (5p) ^^3/2 state 
by the Ba+ cooling laser, the 493 nm and 780 nm laser beams were alternately blocked using a mechanical chopper 
operating at a frequency of 1000 Hz. 

The populations of the (5s) ^5'i/2 and {bp) ^^3/2 levels of the ^'^Rb atoms in the MOT were determined from 
the fluorescence intensity of the atom cloud using the steady-state solution of the optical Bloch equations (OBE) 
for a two level system ^5 . The steady-state populations of the (6s) ^5'i/2, (6p) ^^3/2, and {bd) ^ 1^3/2 levels of the 
Coulomb-crystallised ^^^Ba+ ions were calculated using an 8-level optical Bloch equation model including the effects of 
magnetic fields [271 EH] • This method has been tested and calibrated previously by comparisons with the fluorescence 
yield of single Ca+ ions [23] . 

Because of the large depth of the ion trap (~ 4 eV), the ionic reaction products remained trapped and were 
sympathetically cooled into the Coulomb crystal by their interaction with the remaining laser-cooled Ba+ ions. The 
masses of the different ion species were determined by resonance-excitation mass spectrometry (REMS) as described 
in Refs. [14l[T7]. Molecular-dynamics (MD) methods were used to simulate the laser-cooling fluorescence images of 
the ions in order to determine the composition and kinetic-energy characteristics of the Coulomb crystals [H [29] . 

The BaRb+ complex has the same valence electronic structure as alkali dimers with two external electrons moving 
in the field of the Rb+ and Ba^+ ionic cores. Thus, the electronic potential curves (PECs) and transition dipole 
moments (TDMs) of the BaRb+ molecule can be conveniently calculated using a full configuration interaction approach 
performed in a configuration space generated by the large Gaussian basis sets reported in Refs. j3QH32] . The two- 
electron Hamiltonian involves an effective core potential for the Ba^+ 23 and for the Rb+ [30 ions together with 
core-polarization potentials [35^ ^36 to account for the correlation between core and valence electrons. The spin-orbit 
correction to the PECs has not been considered in the present work. 

The cross sections for RA and radiative charge transfer (RCT) are computed using the same method as in our 
previous study on Ca++Rb reactive collisions [23 . The electronic TDM matrix elements between the initial state 
(a ground-state Ba+ ion colliding with a ground-state Rb atom) and the final states (either a ground-state BaRb+ 
molecule in the case of RA or a pair of free Ba and Rb+ particles for RCT) are integrated over the internuclear 
distance R using the Milne phase amplitude method [371 [SH] . Because of the large mass of the collisional complex, a 
large number of partial waves contribute to the collisions even at the low energies achieved in the present study, so 
that this approach allows for an efficient localization of orbiting resonances (see below). 
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FIG. 2: (a) Demonstration of Ba^ isotope mixing in the Coulomb crystals: (i) Uni-directional laser cooling results in radiation 
pressure on the ^^^Ba^ ions and a localization of the non-fluorescing Ba^ isotopes (predominantly ^^^Ba^, ^"^^Ba^, ^"^^Ba^ 
and ^'^'^Ba^) on the right-hand side of the crystal, (ii) The radiation pressure is cancelled in a bi-directional laser-cooling 
configuration leading to a complete mixing of the isotopes in the crystal, (b) Images of a crystal of Ba^ ions over the course of 
a typical reaction with Rb atoms. Solid white ellipses highlight the approximative boundaries of the crystal of fluorescing Ba^ 
ions, and dashed white ellipses approximately delimit the dark central core of product Rb^ ions, (c) Plot of ln(Vg^+ ^/Vg^+ q) 
against reaction time following pseudo-first-order kinetics. The red line represents a linear regression to the data. See text for 
details, (d) Effective rate constant A; as a function of the population in the excited Rb(5p) ^P^ii state. 



III. RESULTS AND DISCUSSION 



A. Reaction channels and channel-specific rate constants 

Fig. [2] (a) (i) shows a typical Coulomb crystal of laser-cooled Ba+ ions used in the present experiments. Only the 
-^^^Ba+ isotope (72% natural abundance) is laser-cooled. The other isotopes (predominantly ^^''Ba+ (11%), ^^^Ba+ 
(8%), ^^^Ba+ (7%) and ^^^Ba+(2%)) are sympathetically cooled and appear as a dark non-fluorescing region on the 
right-hand side of the crystal. The ^^^Ba+ ions are localized on the left-hand side because of radiation-pressure forces 
exerted by the cooling laser beam propagating along the trap symmetry axis from right to left in the image. When 
the cooling laser is retro-reflected after leaving the trap, the radiation pressure is cancelled and a complete mixing of 
all Ba+ isotopes in the crystal can be observed (Fig. ^ (a) (ii)). 

As reactions of the Ba+ ions with Rb proceeded, the formation of a dark core about the central axis of the 
crystals could be observed which consists of sympathetically-cooled product ions (Fig. |2|b)(ii-iv)). In principle, all 
Ba+ isotopes can participate in the reactions. However, the laser-cooled ^^^Ba+ species reacts preferentially as a 
consequence of the electronic excitation (see below). 

At the low number densities of Rb atoms in our MOT (nR\) ^ 1 x 10^ cm~^), ternary reactive processes which are 
observed in hybrid-trap experiments using denser atom clouds [18 are precluded and the reaction kinetics follow a 
second-order rate law. The constant replenishment of the Rb atoms in the MOT from background vapour results in 
pseudo-first-order kinetics defined by an integrated rate law of the form \^{yBa-^ .tl^Ba-^ ,t^^ — ~kpfot [14l|39j. kpfo 
denotes the pseudo- first order rate constant, t is the reaction time, and VBa+,t=o ^Ba+,t ^ire the volumes occupied 
by the reactant ions at t = and respectively. The second order rate constant k is obtained via k = kpfo/riRi). 

Reaction rate constants were determined by monitoring the decrease of the volume occupied by Ba+ ions as a 
function of the reaction time. For this purpose, the volume V^ot of the entire crystal (i.e., reactant and product ions) 
was estimated by fitting an ellipsoid around its perimeter and subtracting the volume Kore of the central core which 
consists only of sympathetically-cooled product ions (see Figs.|2];b)(ii-iv)). The core of sympathetically-cooled ions 
exhibits in general a more complex shape [40 . However, in order to enable a quick and efficient analysis of the large 
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FIG. 3: (a) Experimental false-colour fluorescence images of Ba^ Coulomb crystals (upper panels) and their molecular- 
dynamics (MD) simulations (lower panels) (i) before and (ii) after reaction with ultracold Rb atoms. Sympathetically-cooled 
ions have been made visible in the MD simulations for clarity. Colour code: blue: ^^^Ba^, green: lighter Ba^ isotopes, yellow: 
^^Rb, red: BaRb^ molecular ions, (b) Resonance excitation mass spectra (REMS) of crystals (i) before and (ii) after reaction. 
See text for details. 



number of images accumulated in the experiments, Vcore was approximated by fitting another ellipsoid to the edges 
of the dark core as indicated in Figs. |2|b)(ii-iv). A plot of ln(yea+,t/^Ba+,^=o) against t is presented in Fig. ^c). 
The data almost perfectly conform to a pseudo- first-order behaviour justifying our approximations for the calculation 
of the crystal volumes. From a linear regression of the data, k was determined to be 2.0(5) x 10~^^ cm^ s~^ for the 
experiment shown. 

By inspection of the energy level diagram Fig. [ijb), it can be expected that the following channels contribute to 
the observed reactions: (i) Ba+(6s)+Rb(5s), (ii) Ba+(5d)+Rb(5s), (iii) Ba+(6p)+Rb(5s) and (iv) Ba+(6s)+Rb(5p). 
Because of the alternate blocking of the Ba+ and Rb cooling laser beams (see sectionM, reactions with simultaneously 
excited collision partners were precluded. The rate constants for the processes (i)-(ivjare designated kg^kd^kp^ and k*. 
The experimentally observed rate constant /c is a weighted average over the contributions from all possible channels 
and can be formulated as [14^: 

^ = ^[iPS,Ca+ ^PS,Rb)ks ^ Pp,Ca+K ^ PD,Ca+kd ^ PP^RbK]^ (1) 

where pij are the steady state populations of the electronic state i of the reaction partner j. We observe a strong 
dependence of k on the populations in the Ba+(6j9) and Rb(5p) states (see Fig. [2](d)). To quantify these dependencies, 
k was measured as a function of the level populations for a set of 19 different laser intensities and detunings from which 
the channel-specific rate constants ki were determined in a multi-dimensional least-squares fit of the data to Eq. Q. 
Some ki were found to be consistent with zero in the fit so that their value was associated with an upper bound given 
by their la uncertainty. The fit yielded /Cg < 5 x 10~^^ cm^ s~^, kp = 2(1) x 10~^^ cm^ s~^, kd < I x 10~^^ cm^ s~^, 
kg = 7(2) X 10~^^ cm^ s~^. The numbers in parentheses denote the la statistical uncertainty of the fitted values. 
The rate constants are affected by an additional systematic error of ~ 50% resulting from the uncertainty in the 
determination of the Rb number density in the MOT. 



B. Reaction products 

Fig. [3] (a) shows images of a representative Coulomb crystal (i) before and (ii) after reaction together with their 
MD simulations. In these experiments, a uni-directional laser-cooling configuration of the Ba+ ions was adopted to 
ensure a better localization of the ions facilitating the analysis of the crystal composition. 

Fig.jsj^b) shows representative resonant-excitation mass spectra of crystals recorded (i) before and (ii) after reaction. 
Whereas only one feature corresponding to the Ba+ reactant ions was observed before reaction, the spectrum after 
reaction shows three peaks assigned to Rb+, Ba+, and BaRb+ motional resonances. The dashed lines in Fig. [s] (b) 
indicate the theoretical single-ion radial excitation frequencies [5 of the different ion species. 

In the MD simulations in Fig. [sj^a), the dark Ba+ isotopes as well as the product ions have been made visible to 
indicate their spatial distribution. The best reproduction of the experimental fluorescence image in (ii) was obtained 
by assuming a composition of 170 ^^^Ba+ (blue), 116 lighter isotopes of Ba+ in their relative natural abundances 
(green), 90 ^^Rb+ (yellow), and 40 BaRb+ (red). From the relative product ion yields, the BaRb+ : Rb+ branching 
ratio (representing an average over all reaction channels) was estimated to be ~1 : 2.25. 

In our previous study on light-assisted Ca++Rb collisions OES], the formation of Rb2" ions was observed. Their 
occurrence was attributed to consecutive reactions of CaRb+ with Rb. Another mechanism leading to their generation 
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could be RA of the Rb+ product ions with Kb atoms from the MOT. In the present case, however, no clear feature 
assignable to Rb2" (which would be located in between the broad BaRb+ and Ba+ resonances in Fig.|(b) (ii)) was 
observed in the mass spectra. We therefore conclude that the formation of Rb2" constituted at best only a minor 
process in the present experiments and was therefore neglected in the modeling of the Coulomb-crystal composition 
after reaction. 



C. Potential energy curves and reaction mechanisms 

Fig. [4] shows the computed non-relativistic potential energy curves of the lowest electronic states of BaRb+ up to 
the Ba"(65)+Rb(5p) asymptote. The PECs for the lowest electronic states exhibit an overall qualitative agreement 
with those reported in Refs. [20l|4T]. In particular, a double-minimum structure for the A PEC which arises 
from an avoided crossing with the C state has also been found in the calculations of Refs. [20l|4T]. In the results 
of Knecht et al. [41 , the top of the barrier appears to be located slightly above the dissociation limit of this state 
both in their non-relativistic and relativist ic calculations, whereas in the current work the barrier height is predicted 
to be slightly smaller than the dissociation energy. Krych et al. [20 demonstrated that the height of the barrier 
dramatically depends on the size of the configuration space used for the calculations. As discussed below, this feature 
has a crucial influence on the RA rate in the lowest reaction channel. 

The entrance channels accessible in the present experiments are indicated in bold. The lowest channel 
Ba+(65)+Rb(55) correlates with the A and a excited electronic states of the BaRb+ molecule. From 
the A -^11+ state, dip ole- allowed radiative processes connecting to the X ^11+ ground electronic state lead to RCT 
forming Rb+ ions and Ba atoms and RA of BaRb+ molecular ions. 

Non-adiabatic transitions induced by the spin-orbit interaction around the curve crossings with the 6, B and c states 
visible in Fig. [4] could in principle lead to non-radiative charge transfer (NRCT). These states, however, correlate 
with Ba(6s5(i)+Rb+ asymptotes which are energetically not accessible at the low collision energies in the present 
experiments so that the NRCT channels are closed. This situation contrasts with the dynamics of the lowest reaction 
channel in the Ca++Rb system in which the energetic order of the analoguous asymptotes is reversed and NRCT was 
found to be the dominant process [Mj |23l |42] . 

The low reaction rate constant kd < I x 10~^^ cm^ s~^ observed for the Ba+(5(i)+Rb(55) can qualitatively be 
explained by the repulsive character of the relevant PECs. The Ba+(6s)+Rb(5p) channel exhibits the largest rate 
constant among the channels observed (/cj = 7(2) x 10~^^ cm^ The fast rate in this channel could be caused by 

radiative processes connecting to the manifold of lower-lying electronic states. Moreover, the collision rate is expected 
to be about 4-5 times larger than in the other channels because of the strongly attractive interaction between the 
charge of Ba+ and the permanent quadrupole moment of Rb in the state pT] . 

We note that the mechanism for the suppression of reaction channels involving excited neutral atoms discussed in 
Ref. [15] does not seem to play an important role under the conditions of the present study. In the present case, the 
Condon point for excitation to the Ba+(65)+Rb(5p) curve is located at an internuclear distance R ^ 1000 ao- At 
smaller distances, laser excitation of the atom is suppressed as a consequence of the shift of the electronic resonance 
induced by the interaction with the ion. However, at collision energies £^co1i/^b ^ 1 K typical for the present study, 
the collision partners can approach over distances exceeding 9000 ao during the lifetime of the excited Rb(5p) ^^3/2 
state so that reactions with laser-excited Rb atoms can readily occur. Similar observations have also been made in 
N2"+Rb collisions at energies down to £^co1i/^b = 20 mK [17]. 

Similar to Ca+-|-Rb [Mj [23], a large rate constant {kp = 2(1) x 10~^^ cm^ s~^) was also obtained for the re- 
action channel with the highest electronic excitation of the ion accessible by the present laser-cooling scheme, i.e., 
Ba+(6p)+Rb(55). The PECs of the highly-excited electronic states correlated with the Ba+(6p)+Rb(55) asymptote 
could not be reliably extracted from the present set of calculations because their pattern was found to be too complex 
to follow them adiabatically from large to small internuclear distances. Further investigations of these states are 
currently under progress. However, the electronic structure at high excitation energies is qualitatively similar to the 
one found in other mixed alkaline-earth-ion alkali-atom systems such as CaRb+ [14 . At high excitation energies, 
CaRb+ exhibits a high density of electronic states correlating asymptotically with excited electronic states of neutral 
Ca. The high density of states promotes opportunities for RA and RCT as well as NRCT mediated by a multitude 
of curve crossings. Apart from CaRb+ [14^, similar effects have recently been invoked to explain the dynamics of 
light- assisted processes in the BaCa+ system [15]. The high density of electronic states occurring at higher excitation 
energies is a common feature of mixed alkaline-earth-ion alkali-atom systems [43^ and it appears that the fast kinetics 
which has been observed in electronically excited channels across several systems [M] [15i i24j is a direct consequence 
of the associated radiative and non-adiabatic dynamics. 
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Internuclear distance R (a.u.) 

FIG. 4: (a) Computed non-relativistic potential energy curves (PECs) of the BaRb^ molecule. The lowest states have been 
labeled using standard spectroscopic notation. Collision channels relevant for the present study have been indicated in bold 
face. States with A symmetry are not shown, (b) Computed transition dipole moments (TDMs) between the ground (X) and 
the two lowest states {A and C) as well as the lowest state (B). The horizontal dotted line represents the energy of 
the initial collisional state in the lowest channel, while the vertical dotted line locates the position of the top of the barrier in 
the A state. 



D. Radiative cross sections and collision-energy dependence 



In the present experiments, the collision energy was dominated by the micromotion of the ions, i.e., their fast 
motion driven by the RF fields |l7^. By contrast, the contribution of the Rb atoms was negligible on account 
of their low temperature (< 150 jjK). The micromotion leads to non-Maxwellian ion velocity distributions which 
were characterised using MD simulations [29 . As the micromotion energy depends on the position of the ions in 
the trap [5^, the average collision energies {Ecou/^b) in the experiments were varied by changing the size and shape 
of the Coulomb crystals |23l[29]. Fig. |5|a) shows experimental (upper panels) and simulated (lower panels) images 
of Coulomb crystals with their collision energy distributions given in Fig. [5|b). Fig. [s] (c) shows the experimental 
(channel- averaged) rate constants k for the crystals in Fig. [5](a) plotted against the average collision energy calculated 
from the distributions in Fig. [5](b). The rate constant was found to be nearly invariant with average collision energy 
within the uncertainty limits, with a slight increase observable at the highest energies. A similarly flat dependence of 
k on EcoW was also observed in the Ca+ + Rb system [I4l|23]. 
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FIG. 5: (a) False-colour fluorescence images of Coulomb crystals of Ba^ ions of varying shapes and sizes ((i)-(v)) and 
their MD simulations ((vi)-(x)). (b) Energy distributions for collisions with ultracold Rb atoms of the crystals shown in 
(a), (c) Dependence of the state- averaged rate constant k on the average collision energy (^coii). Error bars denote the 2a 
statistical uncertainty of at least three consecutive measurements, (d) Theoretical cross sections for radiative association (RA) 
and radiative charge transfer (RCT) for the lowest entrance channel Ba^(6s)+Rb(5s). (e) Theoretical rate constant ks for 
the lowest channel Ba^(6s)+Rb(5s) obtained by velocity-averaging the total radiative cross section over the relative velocity 
distributions. See text for details. 



Fig. 5] (d) shows the theoretical cross sections a for RA (black) and RCT (red) in the Ba+(65)+Rb(5s) entrance 
channel. The cross sections for RA from the A to the X state are predicted to be about an order of 
magnitude larger than for RCT on account of the more favourable Franck- Condon factors for transitions to bound 
vibrational levels of the X state arising from the double- well structure of the A PEC. By comparison, the non- 
relativistic PEC associated with the lowest entrance channel in Ca++Rb collisions exhibits no such feature and 
the relative efficiencies of RCT and RA are reversed [23]. The cross sections exhibit pronounced modulations with 
collision energy identified as shape resonances arising from the trapping of the collisional wave function behind the 
centrifugal barrier at well-defined collision energies. The predicted values for the total radiative cross section in the 
energy interval between 1 and 100 mK (a ^ 10~^^ — 10~^^ m~^) are slightly lower than the experimental estimate 
from Ref. 113 {a ^ 10"^^ - 10"^^ m-^). 

In Fig. p] (e), theoretical rate constants ks for the lowest entrance channel Ba+(65)+Rb(55) are shown which 
were obtained by averaging the total radiative (i.e., RA+RCT) cross sections over the velocity distributions of the 
Coulomb crystals. Because of convergence problems arising from the double- minimum structure of the A state PEC, 
it was not possible to calculate radiative cross sections for high-angular-momentum collisions arising at energies 
Ecow/^B > 400 mK. However, calculations in similar systems have shown that the radiative cross sections show a 
smooth and monotonous energy dependence at low energies [23 . Thus, in order to obtain a theoretical estimate of 
the rate constants, the base line of the total cross section curve (black in Fig. p](d)) was extrapolated over the entire 
energy range relevant for the present study (^coU < 40 K). As shown in Ref. ]23 , the narrow shape resonances only 
insignificantly contribute to the integrated cross section and could therefore be neglected. The rate constant ks for 
the lowest channel was then calculated at the average collision energies corresponding to the experimental data points 
in Fig. [5] (c) by velocity- averaging the extrapolated cross sections over the relative velocity distributions derived from 
the collision-energy distributions in Fig. [5](b). As can be seen from Fig. [5](e), the predicted rate constant ks slightly 
increases with energy over the interval studied. The theoretical values are compatible with the experimental upper 
bound /cg < 5 X 10~^^ cm^ s""*^. The resemblance of the profiles of the channel- averaged rate constant k (Fig. [s (c)) 
and the predictions for the lowest channel (Fig. |5](e)) suggests that the reaction rate shows a similar dependence on 
energy across all relevant channels i.e., Ba+(65)+Rb(55) Ba+(6p)+Rb(55) and Ba+(65)+Rb(5p). 
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IV. SUMMARY AND CONCLUSIONS 



Cold reactive collisions between laser-cooled Ba+ ions and Rb atoms show features similar to other mixed ion- 
atom systems such as Ca++Rb [I4l|23], Yb++Ca [22], Ba++Ca [15] and Yb++Rb [24 : a strong dependence of the 
reaction rate on the electronic state of the reaction partners, the importance of radiative processes and the formation 
of molecular ions by RA. Evidence for RA in these systems has first been established in Ref. [14] by mass spectrometry. 
In the present study, the generation of molecular ions was directly observed by their sympathetic cooling into the 
Coulomb crystals (see Fig. [31(a)). 

As already argued in Ref" 14], light-assisted processes play an important and often dominant role in the chemistry 
of ion- atom hybrid systems. In all of the systems studied so far except Yb++Ca [22 , a pronounced acceleration of the 
rate upon electronic excitation of the reaction partners has been observed. We attribute this observation to the high 
density of electronic states occurring in mixed ion-atom systems already at moderate excitation resulting in multiple 
opportunities for RA, RCT and NRCT. Thus, these processes emerge as a common theme in mixed ion-atom systems 
underlining the rich chemistry that can be observed under the cold conditions achieved in ion- atom hybrid traps. 

On a more general level, hybrid trapping technology has paved the way for the study of ion- atom collisions at 
an unprecedented level of detail by enabling the determination of state-specific rate constants down to millikelvin 
collision energies. The exquisite degree of control and sensitivity achieved in these experiments opens up perspectives 
to gain an improved understanding of ion-atom collision processes that reaches beyond standard pictures such as the 
classical Langevin model [17, .44 and accounts for non-adiabatic, radiative and quant um-collisional effects such as 
orbiting resonances (see section [ill D and Refs. [211 [23]). However, further experimental and theoretical studies are 
required to characterize these phenomena in more detail. 
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